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We consider vector and axialvector mesons in the framework of a gauged linear sigma model with 
chiral U{Nf)ii x U{Nf)L symmetry. For Nf = 2, we investigate the behavior of the chiral condensate 
and the meson masses as a function of temperature by solving a system of coupled Dyson-Schwinger 
equations derived via the 2PI formalism in double-bubble approximation. We find that the inclusion 
of vector and axialvector mesons tends to sharpen the chiral transition. Within our approximation 
scheme, the mass of the p meson increases by about 100 MeV towards the chiral transition. 



I. INTRODUCTION 

The fundamental theory of the strong interaction is quantum chromodynamics (QCD). QCD has 
a local SU{3)c gauge symmetry which determines the interaction between matter constituents, the 
quarks, and gauge fields, the gluons. For massless quarks, the quark sector of the QCD Lagrangian 
also has a global U{Nf)fi, x U{Nf)L chiral symmetry, where Nf is the number of quark flavors. The 
U{1)a anomaly induced by instantons [1] breaks this symmetry expHcitly to U{Nf)v x SU{Nf)A x 
Z{Nf)A, where vector and axial vector symmetries are introduced via V — R-\-L, A = R — L. 
The discrete Z{Nf)A symmetry plays no role for the dynamics and will be omitted in the following. 
Nonzero, degenerate quark masses explicitly break SU{Nf)A, such that the remaining symmetry 
is U{Nf)v- Non-degenerate quark masses further break this symmetry to corresponding to 

baryon number conservation. 

At low momenta ^ Aqcd, where Aqcd ~ 200 MeV is the QCD scale parameter, quarks and 
gluons are confined inside hadrons. Therefore, on a typical hadronic length scale ~ ■'^qcd ^ 
the SU(3)c gauge symmetry of QCD is (at best) of minor importance, and the interactions between 
hadrons are predominantly determined by the global C/(A^/)/j x U{Nf)L chiral symmetry of QCD. 
In the QCD vacuum, the axial U{Nf)A part of the latter symmetry is spontaneously broken by a 
non- vanishing expectation value of the quark condensate {qq) 7^ [2]. According to Goldstone's 
theorem, this would lead to N'j Goldstone bosons. However, since the explicit symmetry breaking 
induced by the U{1)a anomaly reduces the axial symmetry to SU{Nf)A, spontaneous symmetry 
breaking of the latter gives rise to only N'j — 1 Goldstone bosons. These Goldstone bosons acquire 
a mass due to the explicit chiral symmetry breaking by nonzero quark masses. 

At temperatures of the order of ^ (qq)'^^^, the thermal excitation energy is large enough to restore 
the SU{Nf)A symmetry of QCD. If instantons are sufficiently screened [3], this could additionally 
lead to a restoration of the explicitly broken U{1)a- For vanishing quark masses, the high- and the 
low-tcmperaturc phases of QCD have different symmetries, and therefore must be separated by a 
phase transition. The order of this chiral phase transition is determined by the global symmetry 
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of the QCD Lagrangian; for U{Nf)ji x U{Nf)L, the transition is of first order if Nf > 2, for 
SU{Nf)ji X SU{Nf)L X U{l)v, the transition can be of second order if Nf < 2 [4]. If the quark 
masses are nonzero, the second-order phase transition becomes cross-over. 

Lattice QCD calculations predict the chiral phase transition to happen at a temperature Tc ~ 
150 — 190 MeV [5, 6] for zero quark chemical potential /z. The phase transition temperature is 
expected to decrease when fj, increases and vanishes at some value Hc corresponding to quark number 
densities of the order of a few times nuclear matter density. Therefore, the chiral transition of QCD 
is the only phase transition in a theory of one of the fundamental forces of nature which can be 
studied under laboratory conditions: heavy-ion collision experiments performed at the accelerator 
facilities CERN-SPS, BNL-RHIC and, in the near future, CERN-LHC and GSI-FAIR create matter 
which is sufficiently hot and/or dense, such that chiral symmetry is restored. Indeed, the primary 
goal of these experiments is to find evidence for the restoration of chiral symmetry by the creation 
of the so-called quark-gluon plasma, i.e., the phase of QCD where quarks and gluons are liberated 
from confinement. 

When chiral symmetry is restored, the masses of hadrons with the same quantum numbers ex- 
cept for parity and G-parity, so-called chiral partners, become degenerate. Chiral partners are, for 
instance, the sigma and the pion in the (pseudo-)scalar sector, or the p and the ai in the (axial) 
vector sector. A promising signal for chiral symmetry restoration in heavy-ion collisions is therefore 
to study changes of the spectral properties of hadrons in the hot and dense environment [7 9]. One 
of the prime candidates is the p meson. The p meson decays sufficiently fast (and with - for experi- 
mental purposes - sufficiently large branching ratio) into a pair of dileptons which, due to their small 
(since electromagnetic) cross section, are able to carry information from the hot and dense stages of 
a heavy-ion collision to the detector. The CERES and NA60 experiments at the CERN-SPS have 
found convincing evidence for a modification of the p meson spectral function in Pb-|-Pb and In-|-In 
collisions, respectively [10, 11]. 

It is important to clarify whether the modification of the p meson spectral function observed by 
the CERN-SPS experiments is in any way related to chiral symmetry restoration or is merely due to 
many-body interactions in the hot and dense medium. This question can be decided by calculating 
the dilepton production rate from QCD and then using this rate in a dynamical model for heavy-ion 
collisions in order to compute the dilepton spectrum. The low-invariant mass region of the dilepton 
spectrum is dominated by the decay of hadronic states. Therefore, for the calculation of the dilepton 
rate it is more convenient to apply a low-energy effective theory for QCD, featuring hadronic states 
as elementary degrees of freedom and respecting the chiral symmetries of QCD, rather than using 
QCD itself. Since we are interested in the restoration of chiral symmetry at nonzero temperature, 
the low-energy effective theory of choice is a linear sigma model which treats hadrons and their 
chiral partners on the same footing. 

Linear sigma models have been used for quite some time in order to study chiral symmetry 
restoration in hot and dense strongly interacting matter. For instance, Pisarski and Wilczek have 
applied renormalization group arguments to a U{Nf)v x U{Nf)A symmetric linear sigma model with 
scalar degrees of freedom and have drawn important qualitative conclusions regarding the order of 
the chiral phase transition for different numbers of quark flavors [4]. The calculation of hadronic 
properties at nonzero temperature (and density) faces serious technical difliculties. For instance, for 
the following reasons it is impossible to apply standard perturbation theory. First, it turns out that 
the coupling constants of effective low-energy models of QCD are of the order of one, rendering a 
perturbative series in the coupling constant unreliable. Second, nonzero temperature (or density) 
introduces an additional scale which invalidates the usual power counting in terms of the coupling 
constant [12]. A consistent calculation to a given order in the coupling constant then may require a 
resummation of whole classes of diagrams [13]. 

A convenient technique to perform such a resummation and thus arrive at a particular many-body 
approximation scheme is the so-called two-particle irreducible (2PI) or Cornwall-Jackiw-Tomboulis 
(CJT) formalism [14], which is a relativistic generalisation of the $ functional formalism [15, 16]. 
The 2PI formalism extends the concept of the generating functional T[<j)\ for one-particle irreducible 
(IPI) Green's functions to that of one for 2PI Green's functions T[(p,G], where </> and G are the 
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expectation values of the one- and two-point functions. The central quantity in this formalism is the 
sum of all 2PI vacuum diagrams, T2[4>, G]. Any many-body approximation scheme can be derived 
as a particular truncation of r2 [(f), G] . 

An advantage of the 2PI formalism is that it avoids double counting and fulfills detailed balance 
relations and thus is thermodynamically consistent. Another advantage is that the Noether currents 
are conserved for an arbitrary truncation of r2 as long as the one- and two-point functions transform 
as rank-1 and -2 tensors. A disadvantage is that Ward-Takahashi identities for higher-order vertex 
functions are no longer fulfilled [17]. As a consequence, Goldstone's theorem is violated [18, 19]. 
Another consequence is that sum rules of the Weinberg type at zero and nonzero temperature [20] 
are not necessarily fulfilled. A strategy to restore Goldstone's theorem is to perform a so-called 
"external" resummation of random-phase-approximation diagrams with internal lines given by the 
full propagators of the underlying approximation in the 2PI formalism [17]. In this work, however, 
this problem is less severe since we shall only focus on the case of explicit chiral symmetry breaking 
by (small) non-vanishing quark masses. 

Effective theories with chiral U{Nf)ji x U{Nf)L and 0{N) symmetry have already been studied 
within the 2PI formalism in the so-called double-bubble approximation which gives rise to Hartree- 
Fock-type self-energies for the individual particles [18, 19, 21-26]. The 0{N) model has also been 
studied within the 2PI formalism beyond the double-bubble approximation including sunset-type 
diagrams [27]. The chiral models in these works, however, only include scalar and pseudoscalar 
particles. Here, we extend these investigations to vector and axial vector mesons. 

Our ultimate goal is a calculation of the dilepton rate for given temperature (and density). In 
this work, we perform the first step in this direction studying chiral symmetry restoration in the 
mesonic mass spectrum at nonzero temperature. We apply the so-called gauged linear sigma model 
as outlined in Ref. [28]. This model has been used by Pisarski [7] who varied the tree-level mass 
parameters in order to derive qualitative statements about the behavior of meson masses as a function 
of temperature. Here, we extend these studies to a fuU-fiedged self-consistent many-body calculation 
of meson masses at nonzero temperature using the 2PI formalism in double-bubble approximation. 

This paper is organized as follows. The following section is dedicated to a discussion of the 
gauged linear sigma model with U{Nf)ii x U^Nf)^ chiral symmetry, with Nf being the number of 
quark flavors. In Sec. HI we consider the case Nf = 2 in more detail and rewrite the Lagrangian 
explicitly in terms of scalar, pseudoscalar, vector, and axial vector degrees of freedom. We also 
discuss issues related to the breaking of chiral symmetry such as the mixing of pseudoscalar and 
axial vector mesons. In Sec. IV we use this model to calculate the behavior of the meson masses 
and the condensate as a function of temperature. Section V concludes this paper with a summary 
of our results and an outlook for further studies. 

We use the imaginary-time formalism to compute quantities at nonzero temperature. Our notation 

is 



In this section we present the gauged hnear sigma model with U{Nf)ji x U{Nf)L symmetry. The 
Lagrangian reads [28] 




(1) 



We use units 7i = c = fcs = 1. The metric tensor is g^" = diag(l, —1, —1, — 1). 



II. THE GAUGED LINEAR SIGMA MODEL 




(2) 
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Here, the field $ stands for scalar and pseudoscalar hadronic degrees. These are chosen to be 
organized into representations [N^,Nf] under U{Nf)ii x U{Nf)L transformations, such that 

$^4$C/l, (3) 

where J7r,l are unitary matrices acting on the fundamental representation of U(N f)^^^. This allows 
for an explicit representation of $ in terms of a complex Nf x Nf matrix, 

$ = (aa + iTTa) Ta , (4) 

where Ta, a = 0, . . . , — 1, are the generators of U{Nf) in the fundamental representation, and 
TTa, parametrize scalar and pseudoscalar fields, respectively. 
In order to introduce vector and axialvector degrees of freedom, one defines right- and left-handed 
vector fields, 

= {V,^ + A^)T, , L'' = (V7 - A^)T„ , (5) 

where V^, are vector and axialvector fields, respectively. These are treated as massive Yang-Mills 
fields, with field-strength tensor 

BF-^ = d^'E!' - d^R^" - ig [ii", R""] , L^"" = d^'L'' - d'^L^ - ig [L^', L""] . (6) 

Right- and left-handed fields transform as gauge fields under U{Nf)R^L, 

i?" ^ [/^ ( i?^ + - on Ur, ^ ui ( + - dA Ul , (7) 
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while the field-strength tensors transform covariantly under U{Nf)ii^L, 

Rf"" uIr^'^Ur , L^"" uIl^'^Ul . (8) 

The covariant derivative in Eq. (2) couples scalar and pseudoscalar degrees of freedom to right- and 
left-handed vector fields, 

D^^ = d^^ + ig{^L^-R^^) . (9) 

The first line in Eq. (2) contains terms which are invariant under local U{Nf)R x U{Nf)i^ 
transformations. The vector meson mass term renders this symmetry a global one. The de- 
terminant terms represent the U{1)a anomaly of QCD and break the symmetry explicitly to 
SU{Nf)ii X SU{Nf)L X U{l)v, where U{l)v stands for baryon number conservation. The last 
term in Eq. (2) corresponds to the quark mass term in the QCD Lagrangian. Since this term is 

flavor-diagonal, we have H = X^jJ^^ hi2_iTi2_i. For degenerate nonzero quark masses, ho ^ 0, while 
all other /ij2_i vanish. In this case, U{Nf)ji x U{Nf)L is explicitly broken to U{Nf)v- For non- 
degenerate quark masses, also the other /ij2_i have to be nonzero. Exact SU{2)v isospin symmetry 
requires /13 = 0. 

Local U{Nf)ii x U{Nf) l invariance implies universality of the coupling, i.e., the coupling g between 
right- and left-handed vector fields to scalar and pseudoscalar fields is the same as the coupling of the 
vector fields among themselves. Note that the Lagrangian (2) would also be locally invariant under 
U{Nf)ii X U{Nf)L, had we introduced separate coupling constants gR, g^ for right- and left-handed 
vector fields and modified the covariant derivative, -D^$ = d^^ + [{gL^L/^ — 5^!?^$). However, if 
9r 7^ 9l, one obtains non- vanishing parity- violating terms ~ (/ij — 51^ in the Lagrangian, which must 
not occur in viable theories of the strong interaction. 

The assumption of local invariance further restricts possible couplings between scalar and vector 
mesons. Under global transformations, terms like and R^^ would transform as $ itself. Taking 
discrete symmetries into account possible coupling terms would be Tr |$Z/'' — Tr \^L^+R^^\'^ 
and Tr + \R''<P\^) [7]. 
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According to Noether's theorem, continuous symmetries lead to conserved currents. Explicit 
symmetry breaking induces source terms in the conservation laws. For global symmetries, there 
is a simple way to derive the currents and the conservation laws [29]. Consider the symmetry 
transformation to be of the form U = cxp{—i6aTa). Promoting the global symmetry to a local one, 
&a{X), and computing the variation 6£ of the Lagrangian, one can then read off the Noether 
currents and the conservation laws from the identity 

sc = -d^j,^ea-j^d^,ea. (lo) 

For QCD, the U{Nf)ii x U{Nf)L = U{Nf)v x U{Nf)A symmetry leads to the following vector and 
axial-vector currents and the conservation laws: 



dSCqcD 
dSCqcD 



W = -^^^=^m,Ta]q, (11c) 



W = -^^^=ia{ra,M}75g + <5„o|^e,.„^GrGf , (lid) 



where M is the quark mass matrix and G^'' the gluon field-strength tensor. The last term on the 
right-hand side of the conservation law for the axial current cannot be obtained from the variation 
of the classical QCD Lagrangian. It represents the U{1)a anomaly and arises from instantons [1]. 

In the gauged linear sigma model, the vector and axial-vector currents and the corresponding 
conservation laws can be obtained analogously from the variation of the Lagrangian (2) under local 
U{Nf)v X U{Nf)A transformations. The result is 



d6£ ™lTra -70 "^1 Aa 



^^^'^v'" = - ^ = -fabcCTbhc , d^JZ = - ^ = dMhc + 4<5„o c Im [det ($)] , (12b) 

where fabcdabc are the totally antisymmetric and symmetric structure constants of SU{N), respec- 
tively. Since the first line in Eq. (2) is locally invariant imdcr U{Nf)v x U{Nf)A, it cannot contribute 
to the currents or the conservation laws. Then, the vector meson mass term which violates local 
U{Nf)v X U{Nf)A invariance gives rise to the celebrated current-field proportionality (12a) [28]. 

One may wondc;r why the vector current, J^y^, and axial vector current. J^fi- '^rising from scalar 
and pseudoscalar particles do not appear in the expressions (12a) for the vector and axial vector 
currents. They read explicitly 

Jv^. = -J fabc [€(^M0)c - [D^mc] , (13a) 



Ja^. = dabc Wb{D^<l>), - {D^m,] , (13b) 



where 

{D^<l>), = [d^Sac + g {fabcV^ - i46c4)] '^c (14) 

is the covariant derivative. However, employing the equations of motion, these currents are recog- 
nized as simply being parts of the total vector and axial vector currents (12a), 

Ja^ = ^A,--g - fabc {V^Al^ + AiV^^) , (15b) 
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where 

V;, = d^V^-d.V^+gUo{V'^V^ + AlAl) , (16a) 
Al, = d^At-d.A;+gU4V^At + AlV,') (16b) 

are the field strength tensors for vector and axial vector fields, respectively. 

In the following, we briefly discuss consequences of Eqs. (12b). In the chiral limit, ha = 0, and the 
isosinglet vector current Jy^ and the vector currents J7y^, i ^ 0, as well as the axial vector currents 
JTJ^, are exactly conserved. The isosinglet axial vector current receives a contribution from the 
U{1)a anomaly. In the case of explicit chiral symmetry breaking, we have ha ^ 0. However, since 
/ofec = 0, the isosinglet vector current is still exactly conserved. If only /iq 0, but /ij = 0, i ^ 0, for 
the same reason also the other vector currents Jy are exactly conserved. For instance, this is the 
case for = 2 assuming isospin invariancc which requires hi = for i = 1, 2, 3. In contrast, since 
dabo = Sab J the axial vector currents are only partially conserved; the famous partial conservation of 
axial currents (PCAC). 



III. THE CASE Nf ^2 

So far, the discussion of the Lagrangian (2) was valid for an arbitrary number of quark flavours. 
In the following, we restrict ourselves to the case of mass-degenerate up and down quarks. For this 
case, the fields i?^, and can be written in terms of the physical scalar (a. So), pseudoscalar 
[t], tt), as well as vector (w'*, p'^) and axial vector fields (/f , af ) as: 

$ = {aa + ITTa) ta ^ {(J + IT]) to + (do + m) ■ t , (17a) 
i?" = {Va^ + A^J ta = {00^ + /r) to + + «r) • t, (17b) 

L'' = (KT -A':)ta = {CO'' - /r ) to + ip'^ - af) • t . (17c) 
The vector and axial vector mesons enter the Lagrangian through the following terms: 

Tr[(D^$)tD^$] = i[d^a + g{rjh^+Tf-d,^)f 

1 2 

+ 2 i^i^^ - 5(0" htJ. + ao • ai,^)] 
1 2 

+ 2 [^M^O + 9{Pn X ao + 7? Oi^ + TT /i^)] 
1 2 

+ 2 I^/*^ - 5(7? X P/* + o- ai^ + do /i^)] , (18a) 



1 2 

\d^Pv - dvPp. + g{Pp. X + ai^ X dxv)\ 

1 2 
- ^ {dixd\v - d^di^ + g{pfj, X di^ + di^ X Pu)\ , (18b) 

^ Tr[{Ln' + {Rn'] = ^ K + + ff, + O • (18^) 

Note that the w-meson completely decouples from the dynamics. In order to have non-vanishing 
coupling of the u) to the other fields, we would have to include the Wess-Zumino-Witten term [30, 31]. 
Note that, had we included globally U{Nf)v x ?7(A^/)A-invariant terms such as e.g. Tr \^L^^ + R'^^\'^, 
the w-meson would not decouple. 
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Since the w-meson is protected by the U{l)v symmetry, it does not mix with the /i. Consequently, 
we do not expect the lo to become degenerate in mass with the /i when chiral symmetry is restored. 
This is in contrast to the p and the ai, which mix under SU{2)r x SU{2)l transformations: these 
mesons arc expected to become degenerate in mass when chiral symmetry is restored. 

In the phase where chiral symmetry is broken, the a-field assumes a non-vanishing expectation 
value, {(j) = (j) = const. ^ 0. In order to examine the fluctuations around the ground state, we shift 
the (T-field by its vacuum expectation value, a ^ a + (j). After this shift the covariant terms (18a) 
read 

1 2 

+ 2 - 5(0- hij, + (t>fii^ + ao- ai^)] 

1 2 

+ 2 {dfiao + giPt, X ao + ?7ai^ +7r/i^)] 

1 2 
+ -^[di^'^- 9{^^ Pn+(^aii^ + 'l>aii^+aofi,j,)] ■ (19) 

The shift of the cr-ficld leads (after an integration by parts) to the bilinear terms S^/i* and 
5077 • duS^. Physically, these terms correspond to a mixing between the longitudinal component 
of the axial vector mesons and the Goldstone modes arising from chiral symmetry breaking. If the 

mesons were true gauge fields, i.e., mi ~ 0, and if c = iJ = 0, the theory would be invariant under 
U{Nf)ii X U{Nf)A gauge transformations. In this case, the mixing could be removed by an 't Hooft 
gauge-fixing term 

Cgf = -^(5/.^^)' - ^{d^P' f - ^(^M/r + ^9<l>ir- j^id^ar + ^ 5</'vr)^ (20) 

In unitary gauge, ^ ^ oo, the pseudoscalar particles (the Goldstone modes from spontaneously 
breaking U{Nf)v x U{Nf)A to U{Nf)v) would become infinitely heavy, i.e, they are no longer 
dynamical degrees of freedom. 

Since mi ^ 0, the 't Hooft gauge-fixing procedure is not at our disposal. We follow a method 
commonly used in the literature [7, 28, 32], which is a redefinition of the axial fields, 

K = fi^+wd^jj, (21a) 
af = a['' + wd''n, (21b) 

where the new fields ared denoted by a prime (which is later omitted) and w is defined such that 
the bilinear terms mixing the axial vector and pseudoscalar fields are eliminated. 



^ = ^2 ^ ,„.,2 • (22) 



mf -I- {g(j)) 

One could also perform a redefinition of the axial vector fields using covariant derivatives [28, 32], 

/r = + ^ {d^v-gSr-So) , (23a) 
af = ai^'+w [d^n - 5 ( tt x + 5 /f ao)] . (23b) 

Note, however, that the old fields appear in the covariant derivates and couple the set of equations. 
Solving for the old fields in terms of the new one obtains 

/r = {fi'^ + ^" - 9^0 ■ K + (d^^ + X ^)] }) (1 - 9^w^4y' , (24a) 
= {a[^ + w {S'^TT -g[7fxp^ + ao {f"^ + w a^r?)] }) • (l - g^w^SoSo)'' . (24b) 

Note that the term aoOo on the right-hand side in Eq. (24b) is a dyadic product, and consequently the 
corresponding term in parentheses a matrix in isospin space, which has to be multiplied from the left 
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with the preceding isospin vector. The ensuing set of equations couples different isospin components 
of the old and new Oi fields. All this can be avoided by simply performing the redefinition with the 
partial derivative, Eq. (21). It should be noted that such complications did not arise in previous 
treatments [28, 32], because they did not consider the ??, oq, and /i mesons. In that case, the 
redefinition of the ai field simply reads af = a'-^ + w {d'^n + g x if) .It has been noted in Ref. 
[32] that this gives the correct seagull term when coupling the vector mesons to the photon. However, 
in this work we are not concerned with coupling the vector mesons to leptonic currents, and thus 
we restrict the consideration to the simpler version (21) of the redefinition. A different method to 
cope with the problem is to define a non-diagonal propagator [33] . 

After the redefinition (21), the kinetic terms of the pseudoscalar mesons acquire a wave- function 
renormalization, 

1 1 



2 (^/'^a) 



where 



ml + {g(j)) 



(25) 



(26) 



In order to have correctly normalized asymptotic states, we have to redefine the pseudoscalar fields, 

TTa > ZlTa- (27) 

After the redefinition of the pseudoscalar particles the Lagrangian reads 
1 



- [d^a + gZ (ry/i^ + wZiqd^iq + tt • ai^ + wZif ■ d^if) 



1 

+ 2 



{d^i-qf + g"^ {(7 flu + wZad^ri + ao • ai^^ + wZHq ■ d^nf 

- 2g {Zd^T] - g(j)f!^ - g(j>wZd^'r]) (cr/i^ + wZad^r] + clq ■ + wZclq ■ d^n) 

1 2 

+ 2 [^Ai«o + 5 (P/i X «o + ZrjSi^ + wZ'^tjO^n + Ztt/i^ + wZ'^ifd^'q)] 
+\ {^i^Tff + g^ {Ztt X Pi^+ aSin + wZad^Tf + oo/i^i + wZaodi^rif 

- 2g {Zd^TT - g(j)a'( - gcjywZd^n) ■ {Zjf x + aai^ + wZad^Tf + aofif^ + wZaod^ri) 

1 2 1 2 

1 2 

-- {di_iPi, - dupfj, -\- g[pii ^ Pv -\- oi/n X a\u + wZ (ai^ x d^if + 5^7? x di„) + v^Z'^d^ir x d„Tf] } 

1 2 

-- {9^aiy - d^a-iu +9[PnX aii^ + ai/n x p^+wZ {p^ x d^n + d^n x p^)]} 



c + 3( Ai + ^ 



m^ + c+ (Ai + 3y )c6' 



2 
Z^ 
2 



mn + C + ( Ai + y ) 0^ 



mo - c + ( Ai + Y 



2 , -2 



-J (ai + ^) [a' + ZW + al + Z^TT^)^ - (ai + ^) [a' + Z\^ + + Z^^^) 



[ado + Z^Tyit) + Z'^a^'n^ — Z^ (oq • tt) — X24>do ■ [ado + Z^ijifj 



(28) 
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where 



V{4>) = 2 



(29) 



is the classical potential energy density. In the derivation of Eq. (28), we have exploited the fact 
that <j) is the minimum of the potential energy density, i.e.. 







c+ Ai + -f 



ho 



(30) 



Prom the Lagrangian (28), one reads off the tree-level masses for the scalar, pseudoscalar, vector, 
and axial vector mesons. 



TO. 



TO„ 



TOq — c + 3 Ai + 



A, 



A, 



A. 



TO,,, = 



TO^+C+ ( Al + y ^ ^2 

Ai + 3y j 02 



2 2 
TOp = TOi , 

,2 _2 



= TO„^ = TOi + (5f(A) 



(31a) 
(31b) 
(31c) 

(31d) 

(31e) 
(31f) 



Note the factor Z'^ in the definition of the (squared) mass for the pseudoscalar particles. We conclude 

that the wave-function rcnormalization factor Z in Eq. (26) is equal to the ratio of the trcc-lcvcl 
masses of the ai and the p meson, Z = ma^/mp. In the framework of this model, the difference in 
mass of the two mesons is due to the Higgs effect. 

Using the expression for the pion mass, we obtain from Eq. (30) that /iq = Z~'^rn^(f). Using the 
renormalization of the pion field (27), the second Eq. (12b) reads df^J^ = Z~^m'^(j)'Ki. From the 
PCAC relation dp,J^^ = /ttTO^ ttj we then conclude that the vacuum expectation value of the scalar 
field is 



(t> = z u = U 



TOq. 

rrin 



Note that the KSFR relation [34] 



TO„ 



(32) 



(33) 



makes a definite prediction for Z, and thus for the relation between the mass of the p and the ai: 
Eqs. (26), (32), and (33) imply that Z = \f2, or that = \f2rap. For vfip ~ 770 MeV, this means 
that ~ 1090 MeV, somewhat smaller than the value quoted by the PDG [35], nia^ ~ 1230 MeV. 

The vector and axial vector mesons are massive spin-one particles, i.e., they have three physical 
degrees of freedom. In order to have an invcrtiblc inverse trcc-lcvcl propagator for these particles, 
we have to promote the unphysical fourth degree of freedom of the vector and axial vector fields to 
a physical degree of freedom by adding a Stueckelberg term, 



The unphysical degrees of freedom can be removed by taking the limit ^ 
calculation. 



(34) 

at the end of the 
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FIG. 1: Two-loop order diagrams in V2. 



IV. THERMAL MESON MASSES 



In this section we investigate chiral symmetry restoration at nonzero temperature in the framework 
of the model (28). We explicitly compute the condensate and the meson masses as functions of 
temperature. As mentioned in the introduction, wc use the 2PI rcsummation scheme, or ^-functional 
[16], or CJT formalism [14]. The generating functional for 2PI Green's functions is 

r[cl>,a\ = S[ct>] + ^TrlnG-^ + ^Tr{D-^G-l)+T2[^,G] , (35) 

where </> and G are the expectation values of the one- and two-point functions in the presence of 
external sources, S [(/)] is the tree-level action, is the tree-level propagator and r2 [(f), G] iw the 
sum of all 2PI vacuum diagrams with internal lines given by G. If the system is translationally 
invariant, it suffices to consider the effective potential V = —TT/fl, where T is the temperature and 

n the 3- volume of the system. 

For the Lagrangian (28) the cfFoctive potential is given by 

V [<^, Gi] = V{<P) +Ia + Ir, + Slao + 3/. + /a, + //, + 3/^ + 34, + V2 [(/., Gi] , (36) 

where V2 = -TT2/O., and 



= \ [[\nGr\k)+V-\k)G,{k)-l] 



(37) 



with the inverse tree-level propagators for scalar particles 'D~^ = —k"^ + mf, i = a,ri, qq, tt, and for 
vector particles V^^^ = — (A;^ — mf)g^i, + (1 — Ci^n^v^ i = u, fx, p, Oi. The condensate and the full 
propagators Gj are determined from the stationarity conditions for the effective potential (36), 

The last equation is identical to the Dyson-Schwingcr equation for the full propagators, 

G-^=D-'+lii, (39) 

where the IPI self-energy Ilj is given by 

To lowest order in the loop expansion, i.e., to two- loop order, two types of diagrams contribute to V2. 
The diagram in Fig. 1, made of a single four-point vertex, is called the double- bubble diagram. The 

other one, made of two three-point vertices, is called the sunset diagram. In Eq. (40), the functional 
derivative of V2 with respect to a full propagator is diagrammatically equivalent to cutting a line 
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Q 




FIG. 2: IPI self-energy U. 



corresponding to that propagator. This gives rise to the self-energy diagrams shown in Fig. 2. In 

contrast to the sunset diagram, the double-bubble diagram does not produce an imaginary part in the 
IPl self-energy when analytically continuing the latter to real energies. Thus, it gives only constant 
(temperature-dependent) corrections to the mass and the wave-function renormalization constant 
of the particle. Neglecting the sunset terms is, of course, still a self-consistent approximation in 
the framework of the 2PI formalism. The advantage is that it tremendously reduces the numerical 
eifort: instead of solving a set of coupled complex integral equations for the self-energy functions 
Il{ko,k) for a given temperature, one only has to solve a set of coupled real fixed-point equations 
for the masses and wave-function renormalization constants. 

For the Lagrangian (28), wc explicitly list V2 in double-bubble approximation in Eq. (Al) of App. 
A. As detailed in that appendix, we make a further approximation to simplify the structure of the 
vertices corresponding to interactions of vector and axial vector mesons with themselves and with 
pions. This approximation yields degenerate dispersion relations for transverse and longitudinal 
(axial-) vector degrees of freedom. The resulting self-energies are given in Eqs. (A6), and the set of 
coupled fixed-point equations in Eqs. (A7). 

The parameters of the model are fixed to the vacuum masses and the pion decay constant. Without 
vacuum fluctuations, vacuum and tree-level masses are identical, and we have 



1 



2 




(41a) 



2(^^)2 L' 



-Wan + 




(41b) 



(41c) 



c 



ml -ml 

2^2 



(41d) 




(41e) 
(41f) 




(41g) 



with 



m„ = 400- 1400 MeV 



m,, = 547.75 McV, 

TOao = 985.1 MeV, 

m^ = 138.04 MeV, 

mp = 768.5 MeV, 



ma, = "^-/i = 1230 MeV 
U = 91.9 MeV. 



In Figs. 3 and 4 we show the temperature dependence of the condensate and the scalar meson masses 
without the influence of vector mesons, i.e., setting <? = (and Z = 1). In Fig. 3, we choose a vacuum 
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sigma mass of nia = 441 McV [36]. This leads to a cross-over phase transition at Tc 193 MeV, in 
good agreement with the lattice QCD calculations of Ref. [6]. 

Taking a vacuum sigma mass of nia = 1370 MeV (which is another candidate for a 0++ state [35]), 
the phase transition is of first order at a critical temperature of Tc ~ 238 McV, cf. Fig. 4. Within 
this scenario with only scalar, but no vector mesons, we obtain a second-order phase transition for 
a critical sigma mass of about m* ~ 965 MeV at a temperature of Tc 222 MeV. 

In Figs. 5 and 6, wc show the meson masses, the condensate, and the wave-function rcnormaliza- 
tion constants as functions of temperature including vector meson degrees of freedom. For = 441 
MeV, cf. Fig. 5, we observe a cross-over transition at about the same temperature, Tc 195 MeV, 
as in the pure scalar case. However, the inclusion of vector mesons leads to a more rapid cross-over. 
This is a typical effect within chiral linear a models when including more dynamical degrees of 
freedom, e.g., strange and charmed scalar mesons [22]. Including vector mesons this effect is even 
stronger. 

The increasing steepness of the cross-over transition indicates that the inclusion of vector mesons 
brings one closer to a second-order critical point. This can also be seen comparing the value of 
the a mass at the respective transition temperature in Figs. 3 and 5. In the latter, the a meson is 
lighter by about 100 MeV. At a second-order critical point, the a should become completely massless 
[37], indicating a diverging correlation length. Moving closer to a second-order critical point by the 
inclusion of vector degrees of freedom is advantageous from the following point of view. Let us 
consider the second-order critical endpoint of the line of first-order phase transitions associated with 
the restoration of chiral symmetry, or shortly, the QCD critical point. Usually, in chiral models of 
QCD, this point is located at temperatures which are too small and quark chemical potentials which 
are too large [38] as compared to lattice QCD calculations [39]. Our calculations arc performed 
at /X = 0. The fact that the inclusion of vector mesons brings one closer to a second-order phase 
transition indicates that the QCD critical point moves closer to the T axis in the T — /x plane. We 
thus conjecture that the inclusion of vector mesons in linear a models improves the agreement with 
lattice QCD. 

Note that the p mass increases, and the oi mass decreases, as a function of temperature, until 

these mesons become degenerate in mass beyond the chiral symmetry restoring transition. This is a 
necessary requirement for chiral symmetry restoration, as argued in Ref. [7]. The question remains, 
whether the mass of the p increases to meet that of the ai, or whether both decrease (or, for that 
matter, increase) with T, before becoming degenerate. In our model, the increase of Mp with T 
is not surprising, because it consists of two separate contributions, the tree-level mass parameter 
nil, and the contributions from thermal (tadpole) fluctuations. While the former is constant as a 
function of temperature, one can show that the latter are always non-negative, therefore the mass 
of the p has no choice but to increase with temperature. This would be different, though, had we 
chosen to generate the tree- level mass dynamically through chiral symmetry breaking [32]. In this 
case, the tree-level contribution would decrease as the condensate melts, while the thermal tadpole 
contributions would increase. (To our knowledge, this argument was first given in Ref. [8].) As 
the condensate goes to zero in the process of chiral symmetry restoration, only the latter survive 
and lead to an increase of Mp with T. Nevertheless, in this scenario, at intermediate temperatures 
one cannot exclude that the p meson mass first decreases with T, assumes a minimum, and then 
increases again. The minimum value of Mp may well be zero, as advertised by models with hidden 
gauge symmetry [40]. 

The pion and 77 wave-function renormalization constants are unity at T = 0, and then increase by 
about 10 % (for the pion) or 20 % (for the t]) when approaching the chiral phase transition. They 
sharply drop back to one when the transition temperature is exceeded. 

In Fig. 6 we show the thermal masses, the condensate, and the wave-function renormalization for 
rricr = 1370 MeV. In this case, there is a first-order phase transition at Tc ~ 258 MeV. For the sake 
of completeness we note that the second-order phase transition now occurs for a critical sigma mass 
m* ~ 795 MeV at a temperature of Tc 220 MeV, i.e., at smaller values of and Tc than without 
vector mesons. 
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FIG. 3: Masses for scalar and pseudoscalar mesons (left panel) and chiral condensate (right panel) as a 
function of temperature for the case g = 0, nia = 441 MeV. 




T [MeV] T [MeV] 



FIG. 4: The same as in Fig. 3, for rua = 1370 MeV. 



V. CONCLUSIONS 



In this work, we have investigated an effective theory for QCD with chiral U (2)/^ xU{2)l symmetry. 
Besides scalar and pseudoscalar degrees of freedom, we included vector and axial vector mesons. 
These were coupled minimally to the scalar and pseudoscalar mesons. The ensuing local U (2)i^ x 
U{2)l symmetry is explicitly broken to the corresponding global symmetry by introducing a mass 
term for vector and axial vector mesons. Without introducing a Wess-Zumino-Witten term, it turns 
out that the ui meson completely decouples from the theory. 

In the Goldstone phase, a non-vanishing chiral condensate leads to terms bilinear in the axial 
vector and pseudoscalar fields. Lacking the freedom to choose a convenient (i.e., 't Hooft) gauge, 
in order to get rid of these terms one can perform a redefinition of the axial vector fields with a 
subsequent redefinition of the pseudoscalar degrees of freedom. 

The resulting Lagrangian was employed to derive a set of coupled Dyson- Schwinger equations 
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50 100 150 200 250 30O 350 50 100 150 200 250 300 350 400 
T [MeV] T [MeV] 



FIG. 5: Masses for scalar and pseudoscalar mesons (upper left panel), for vector and axial vector mesons 
(lower left panel), the chiral condensate (upper right panel), and the pion and 77 wave-function renormaliza- 
tion constants as a function of temperature for the full model including vector and axial vector degrees of 
freedom; m,^ = 441 MeV. 

via the 2PI formalism. In this work, we approximate the 2PI effective potential by including only 

the so-called double-bubble diagrams. This leads to constant, temperature-dependent shifts of the 
masses and, in the case of the pseudoscalar mesons, the wave-function renormalization constants. 

We did a control calculation without vector and axial vector meson degrees of freedom in order 
to study the influence of the latter on the chiral phase transition. We found that, without vector 
mesons, and with a vacuum sigma mass nia = 441 MeV, there is a smooth cross-over phase transition 
at a critical temperature of Tc ~ 193 McV. In the case with vector mesons, this cross-over transition 
becomes sharper, however, the transition temperature remains approximately the same, Tc ~ 195 
MeV. 

Increasing the mass of the sigma, the phase transition becomes sharper and turns into a first-order 
transition beyond a critical sigma mass to* . At to* , the transition is second of order. For the case 
without vector mesons, this happens for to* ~ 965 MeV, with a critical temperature of Tc 222 
MeV. For the case with vector mesons, this happens already for a smaller critical sigma mass, 
TO* ~ 795 MeV with a critical temperature of Tc ~ 220 MeV. The fact that to* becomes smaller 
when including vector mesons can be understood from the fact that, for a given rricr, the transition 
was found to become sharper as compared to the case without vector mesons. This indicates that 
the inclusion of vector mesons brings one closer to a second-order critical point. 

This has consequences for the location of the QCD critical point in the T — fi phase diagram. 
Without vector mesons, chiral models of QCD usually predict this point at temperatures which 
are too small and at quark chemical potentials which are too large as compared to lattice QCD 
calculations. Since our above mentioned results arc obtained at /i = 0, the inclusion of vector 
mesons would move this point closer to the T axis and thus would improve the agreement with the 
lattice results. 
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FIG. 6: Same as in Fig. 5; nia = 1370 MeV. 



In our calculations, the mass of the p meson is a monotonously increasing function of tempera- 
ture. This is not surprising; in our model, the mass of the p essentially consists of two contributions, 
the tree-level mass plus thermal fluctuations from tadpole diagrams. The former is always con- 
stant, while the latter are positive definite and, in general, increasing functions of temperature. A 
qualitatively different scenario would be one where the tree-level mass is also generated dynamically 
through spontaneous symmetry breaking. A self-consistent treatment of this case within the Hartree 
approximation could be a topic of future studies. However, the generic feature of chiral symmetry 
restoration, namely that the p and the ai become degenerate in mass at the chiral phase transition, 
should still hold for this case. 

The ultimate goal is to find signatures for chiral symmetry restoration in hot and dense nuclear 
matter as created in heavy-ion collisions. In this context, a promising observable is the invariant 
mass spectrum of dileptons, since they carry direct information from the hot and dense stages 
of the evolution of the system. Of particular interest is the contribution of the p meson to the 
dilepton spectrum, since it has a large (vacuum) decay width and therefore decays inside the medium. 
Changes of the spectral properties of the p meson potentially related to chiral symmetry restoration 
can thus be directly seen in the dilepton invariant mass spectrum. 

With this in mind, the present study should be extended along several lines. In the double-bubble 
approximation to the 2PI effective action, the quasiparticle self-energy has no imaginary part. Thus, 
in this approximation the quasiparticlcs cannot decay. A natural next step is therefore to include 
sunset-type diagrams in the 2P1 effective action, which lead to nonzero imaginary parts for the 
self-energy of the quasiparticles and, in turn, to a nonzero decay width. Another possible project is 
the self-consistent calculation of meson masses for the gauged linear sigma model extended to the 
three-flavor case [28] . From previous studies of dilepton production we know that baryons play an 
important role [9]. Therefore, it is mandatory to include nuclcons into the present model. Of course, 
a linear representation of chiral symmetry requires to also add the chiral partner of the nucleon as an 
explicit degree of freedom [41, 42]. Since the chiral partner of the nucleon is most likely heavier than 
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the A resonance, the latter is more abundant in a hot and dense system. It is therefore unavoidable 
to also include spin 3/2-resonances and their chiral partners [43]. 
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APPENDIX A: EFFECTIVE POTENTIAL AND GAP EQUATIONS 



In this appendix, we simplify our notation by denoting full propagators of a given field by the 
field itself, for instance Ga{k) = <j{k), G'^''{k) = p^''{k) etc. The trace over a vector propagator is 
denoted by the corresponding field without Lorentz indices, for instance p^(A;) = p{k). With this 
notation, the generating functional in double-bubble approximation reads: 
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At this point we decompose the terms in V2 derived from vectorial vertices in the following 
manner. Let Gy" denote a vector or axial- vector propagator (with trace Gy, according to our above 
convention). We then decompose a rank-2 tensor structure B^'^ as follows, 



4 ^ ' 4 ^ 



(A2) 
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and apply this to Gy' as well as k^k'^. We then neglect terms involving SB'^" . This results in the 
replacements 



Gv{k) 



- f G^^(fc) / Gv^Al) 

Jk Jl 



1 2 



Gvik) 
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k^Tr{k) 



(A3c) 



In this way, the vector and axial- vector meson self-energies become proportional to 5^'^, i.e., trans- 
verse and longitudinal self-energies become identical. 

The condensate equation with the approximated vectorial vertices reads 
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(A5a) 
(A5b) 

(A5c) 
(A5d) 
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The approximation (A3) has the distinct feature that the vector-meson self-energies are propor- 
tional to the metric tensor, II^" = g^^'^Ili, i = f\,p,a\. The self-energies with the approximated 
vectorial vertices are 
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Finally, the coupled set of fixed-point equations for the masses and wave-function renormalization 
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constants is given by 
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which has to be solved together with Eq. (A4). A consistent ansatz for the full propagators in 
momentum space is 
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The tadpole integrals emerging after the Matsubara summation are 
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fl,P,ai 



where Lu-i_ = (P + MfY/"^ for i = a,ao, fi, p,ai, while uji = (P + Mf/Zf)^/^ for z = tt, 77, and 

= + A/.f/i;)^''^ for z = /9, oi; nB(w) = [exp(ij/T) - is the Bose-Einstein distribution. For 
the actual computation, we renormalize the tadpole integrals in such a way that divergent vacuum 
contributions vanish. Other renormalization prescriptions (yielding nonzero vacuum contributions) 
are possible, but do not lead to qualitatively different results [19]. Note that, after taking the limit 
^ — > 0, only the three physical vector meson degrees of freedom contribute to the tadpole integrals. 
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